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Abstract. In this paper we are concerned with the convergence rate of so- 
lutions to the three-dimensional turbulent flow equations. By combining the 
U'-L'^ estimates for the linearized equations and an elaborate energy method, 
the convergence rates are obtained in various norms for the solution to the 
equilibrium state in the whole space, when the initial perturbation of the 
equilibrium state is small in //"^-framework. More precisely, the optimal con- 
vergence rates of the solutions and its first order derivatives in L'^-norm are 
obtained when the L*'-norm of the perturbation is bounded for some p G 
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1. Introduction 

We consider in this work the turbulent flow equations for compressible flows on 
pt + div(/9u) = 0, 

{pu)t + div(pu (g) u) - Au - Vdivu + Vp = -|V(pfc), 
{ph)t+AiY{puh)-^h=^^+Sk, 
(pk)t + div(pufc) — Afc = G — pe, 
{pe)t + dw{pue) -Ae=^^^, 
^ (p, u, h, k, e){x, t)\t=o = {po{x),uo{x),ho{x),kQ{x),eo{x)), 

with Sk = Hw- + 7r^)- + ^W^^, G = §^[ile{¥- + 1^) " 

^ ^^^oxj oxi ' 3 oxk^oxj axj oxj ' oxj ^'^^axj oxi ' 

lSij{pk + Mefj^)], where % is given by 6ij = if i 7^ = 1 if z = j, dynamic 

viscosity p and the eddy viscosity /xt are positive constants satisfying fi + fit ^ Me, 
and Ci, C2 are also two adjustable positive constants. 

Here p, u, h, k and s denote the density, velocity, total enthalpy, turbulent 
kinetic energy and rate of viscous dissipation, respectively. The pressure p is a 
smooth function of p. In this paper, without loss of generality, we have renormalized 
some constants to be 1. The system (jl.ip is formed by combining effect of turbulence 
on time-averaged Navier-Stokes equations with the k-e model equations. 

All flows encountered in engineering practice, both simple ones such as two- 
dimensional jets, wakes, pipe flows and flat plate boundary layers and more com- 
plicated three-dimensional ones, become unstable above a certain Reynolds number. 
At low Reynolds numbers flows are laminar. Flows in the laminar regime are de- 
scribed by the continuity and Navier-Stokes equations which have been studied by 
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many people [IllllEllelEllHlElIIOllTlllISlIIHlinillllESl^ 



At high Reynolds numbers flows are observed to become turbulent. A chaotic and 
random state of motion develops in which the velocity and pressure change con- 
tinuously with time within substantial regions of flow. More precisely, at values of 
the Reynolds number above Rccrit a complicated series of evens takes place which 
eventually leads to a radical change of the flow character. In the final state the 
flow behavior is random and chaotic. The motion becomes intrinsically unsteady 
even with constant imposed boundary conditions. The velocity and all other flow 
properties vary in a random and chaotic way. Turbulence stands out as a proto- 
type of multi-scale phenomenon that occurs in nature. It involves wide ranges of 
spatial and temporal scales which makes it very difficult to study analytically and 
prohibitively expensive to simulate computationally. Many, if not most, fiows of 
engineering significance are turbulent, so the turbulent flow regime is not just of the- 
oretical interest. Up to now, although many physicists and mathematicians studied 
turbulent flows, there are not any general theory suitable for them. Fluid engineers 
need access to viable tools capable of representing the effects of turbulence. 

This paper is devoted to study decay rates for the system and proves the 
optimal convergence rates of its solutions under suitable assumptions. Bian-Guo [3] 
has obtained the global existence of smooth solutions to the system (jl.ip under the 
condition that the initial data are close to the equilibrium state in if '^-framework. 
More precisely, this result is expressed in the following. 

Proposition 1.1. Assume that initial data are close enough to the constant state 
(p, 0, 0,fc, 0), i.e. there exists a constant such that if 



then the system admits a unique smooth solution (p, u, h, k, s) such that for 

any t £ [0, oo), 



where C is a positive constant. 

Based on this stability result, the main purpose in this paper is to investigate 
the optimal convergence rates in time to the stationary solution. We remark that 
the convergence rate is an important topic in the study of the fluid dynamics for 
the purpose of the computation [131 IE] ■ The main idea in this paper is to combine 
the U'-L'^ estimates for the linearized equations and an improved energy method 
which includes the estimation on the higher power of L^-norm of solutions. By 
doing this, the optimal convergence rates for the solutions to the nonlinear problem 
()l.ip in various norms can be obtained and are stated in the following theorem. 

Theorem 1.2. Let 5q he the constant defined in Proposition There exist 

constants 5i G (0, (5o) and C > such that the following holds. For any 5 < 5i, if 



\{po - p,uo,ho,ko - fc,eo)||if3(R3) < 6o, 



(1.2) 




(1.3) 



< C'IKpo - p,uo,ho,kn ~ k,eo)\\H3, 



II (Po 



P, UQ,ho, fc( 



^0 



fc,£o)||_fi-3(R3) < S, 



(1.4) 



and for some p G [1, f ); 



p,Uo,ho,ko 



k,eo e LP{R^), 



(1.5) 
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then the smooth solution {p, u, h, k, e) in Proposition li.il enjoys the estimates for 
t e [0, oo), 

\\{p- p,u,h,k-k,e)\\,<C{l + t)-^^P^^'°\ 2<q<6, (1.6) 

\\{p-p,u,h,k-k,e)\\^<C{l+ty''^P'^-^^\ (1.7) 

\\V{p~ p,u,h,k-~k,e)\\H'^ <C{l + t)-^^P^^-^^\ (1.8) 

\\{puUuht,kuet)h<C{l + t)-''^P^^-'^\ (1.9) 
where a{p, q; I) are defined by 

a(p,g;0 = |(---) + ^. (I-IO) 
2 p q 2 

Remark 1.3. (jl.Sp shows that the perturbation of initial data around the constant 
state (p, 0, 0,fc, 0) is bounded in iP-norm, for some p & [1, |), which need not be 
smalL 

Remark 1.4. The Unearized equations of ()1.1|) around the constant state {p, 0, 0, k, 0) 
take the foUowing form: 

at + 7divw = 0, 

Vt + "fVa — XAv — AVdivf = 0, 

ht-XAh^O, (1.11) 
mt — XAm = 0, 
et - XAe = 0, 

where 7, A are positive constants which will be given precisely in Section 2. Com- 
pared to the decay estimates of the solutions to the above linearized equations by 
using Fourier analysis [22] stated in Lemma 2.1 in the next section, Theorem 11.21 
gives the optimal decay rates for the solution in L'-norm, for any 2 < q < 6, and its 
first order estimates in L^-norm. Note that the convergence rates of the derivatives 
of higher order in L^-norm and the solution in L°°-norm are not the same as those 
for linearized equations. 



Remark 1.5. We mainly use the method in [12] to prove Theorem 11.21 But our 
problem is much more difficult because of the strong coupling between velocity, 
total enthalpy, turbulent kinetic energy and rate of viscous dissipation. Moreover, 
our result is better than that in [12| . Here we assume L^'-norm of initial data 
Po — p, Uq, ho, kQ — k, Eq is bounded, for somep e [1, |), instead of \\{pQ—j5, uq, Hq, kQ — 
fc,eo)||Li < 00. 

Notation: Throughout the paper, C stands for a general constant, and may 
change from line to line. The norm ||(A, i?)||jf is equivalent to H^Hjc + and 
||A||xnF = Pllx + ljAlly. The norms in the Sobolev Spaces i?™(E3) and W^™'?(R3) 
are denoted respectively by || • and || • for m > 0, g > 1. In particular, 

for TO = 0, we will simply use || ■ ||p. Moreover, (■, ■) denotes the inner-product in 
L2(R3). Finally, 

V =idi, 82,83), d, = d,^, z = l, 2, 3, 
and for any integer Z > 0, V'/ denotes all derivatives up to Z-ordcr of the function 
/. And for multi-indices a, /3 and .J 

a = (ai,Q;2,Q;3), /3 = (A, /52, /Ss), ■? = (?i,6,6), 
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we use 



and C£ 



a - p\{a-py. when /3 < a. 



1=1 



2. Preliminaries 



We will reformulate the problem (jl.ip as follows. Set 



7 = 

Introducing new variables by 

a = p — p, V 



p'{p) + k, X = -. 

P 



1 



, h ~ h, m 



k, e 



(2.1) 



the initial value problem (jl.ip is reformulated as 

at + 7divw = Fi , 
i;t + 7Va — XAv — AVdivu = _F2, 

ht - XAh = Fg, 

m( — AAm = F4, 
et-\Ae = F5, 

^ (a, w, h, m, s){x, 0) = (oq, wq, /iqi ™0: Eo) —J' as |x| — ?> 00, 

where 

Fi = — 7Adiv(aM), 

i., . ,i - i,(A„ + Vdiv.,) - - m + |!^±|1 - p)Va 

p p ^X^ a + p p 3(a + pj 3p' 

F3 = (- - i)A/i - 7Ap'(a + p)divw H ^-^S'^ - 7AU • V/i, 

p jO a + p 

= (- - i)Am H - e - 7AV • Vm, 

p p a + p 

F5 = (- - -)Ae + _ - - -^Sk - jXv ■ Ve, 
p p (a + p){m + k) 171 + k 

with iS*^ and the corresponding Sk and G in the variables of (a, u, ft., m, e). 
Let [/ = (a,w), C/o = (ao, wo), = (J^i, ^^2), 

7div 
7V AA + AVdiv 

and E{t) be the semigroup generated by the linear operator A, then we can rewrite 
the solution for the first two equations of the nonlinear problem as 

(2.3) 



(2.2) 



A 



U{t) = E{t)Uo + [ E{t~ s)Fds. 
Jo 



The semigroup E(t) has the following properties on the decay in time, which can 
be found in [211 122] and will be applied to the integral formula ()2.3p . 
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Lemma 2.1. Let I > be an integer and l<p<2<q< oo. Then for any t >0, 
it holds that 

\\VE{t)Uo\U<Cil + t)-<P''^-''^\\Uoh.nH'. (2.4) 
with <7{p, q; I) defined by (|1.10[) . 

To treat the last three equations, we introduce the semigroup S{t) generated by 
AA, then (Emilia -(EH) 5 become 



h{t) = S{t)hQ 



f S{t~ s)Fsds, 
Jo 



m{t) = S{t)mo + S{t- s)Fids, (2.5) 



e{t) = Sit)eo + f Sit- s)F5ds. 
Jo 



(2.6) 



We state the large-time behavior of solutions to the last three equations of the 
system (|2.ip as the following lemma which can be obtained by direct calculation or 
can refer to [^ . 

Lemma 2.2. For the solution {h,m,e) of the last three equations of the system 
()2.1|) with Cauchy data h{x,Q) = Hq, m{x,0) = mo, e(x,0) — Eq, there exists a 
constant C such that 

\\{VK V'm, V'e)||, < C(l + t)-(f'9;0 

+ C [ il+t~s)--^P^''-^'^\\{F3,F^,F,)\\p, 1^0, 1, 
Jo 

for any t > 0, ^ < p,q < +oo, as well as a is defined by (|1.10p . 

For later use we list some Sobolev inequalities as follows, cf. pi 111). 

Lemma 2.3. Let 51 C M"^ be the whole space R"^, or half space R+ or the exterior 
domain of a bounded region with smooth boundary. Then 

fi) ll/IU«(0) < C||V/|U.(a), for f e H\n). 

(n) \\f\\L^in)<C\\f\\HHn), for2<p<Q. 

(^) ll/llco(n) < C||/||vvi,p(o) < C||V/||Hi(n), for f e W'i^). 

H \J^f-g-hdx\ < £||V/||2, + ^\\g\\l4h\\l., for e > 0, /, g e H\n), 

h e L^{VL). 

(v) \S^f-g-hdx\< e||.9||i. + f|lV/||^,||/j||2,, for e > f e H\n), g, h G 
L2(f7). 

Finally, the following elementary inequality |12l will also be used. 

Lemma 2.4. If ri > 1 and 7'2 £ [0,ri], then it holds that 
t 

(l + t-sy^il + sy^ds < Ci(ri,r2)(l + i)-''% (2.7) 



where Ci(ri,r2) is defined by 



0^2 + 1 

Ci(ri,r2) = -. (2.8) 

ri - 1 
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3. Basic estimates 

In this section we shall establish two basic inequalities for the proof of the optimal 
convergence rates in section 4. One inequality is the decay rate of the first order 
derivatives, while the other is the energy estimate. 

Lemma 3.1. Let W = {U, h,m^ e) be the solution to the problem (|2.ip . then under 
the assumptions of Theorem \1.2l we have 

\\VWh<CEoil + t)"''^P'^'^^ +C f\l + t~s)-''^>''''''^\\VW\\H2ds, (3.1) 

Jo 

with Eq = \\W{0)\\LPnH^ = ||(C/o,/io,TOo,eo)||LPnffi = \\{ao,vo, ho,mo, Eq^lp^h^ 
and 1 < p < I . 

Proof. Let / = 1 in (|2^ . we then have from ((23)) 

\\VUm2<CEo{l + t)-'^^P'^-''^ +C f\l + t-s)-'^^P^'^''^\\iF,,F2)h.nmds, 



which together with (|2.6p implies that 

\\WW{t)h<CEo{l + t)-'^^P^'-^^nC f\l + t-s)-^^P^^-^^^[\m,F2)hpnH^ 

Jo 

+ ||(F3,F4,F5)l|p]ds. 

(3.2) 

For ||(Fi,F2)|liPnffii estimate as follows. For 1 < p < |, the term in Fi can be 
estimated as 

II - 7Adiv(a«)||j, < Cmav'Wp + \\ady\\p) 

<C||f||^||Va||2 + C||Vz;||2||a||^ 

2—p 2— p 

<C<5||(Va,ViOI|2. 
With the help of Holder inequality and Lemma 12.31 it holds that 

II -7Adiv(aiO||Hi < Cmav'h + Wad.v'h + || Vdiv(a«)||2) 
<CS\\{\/a,\/v)\\Hi. 
Similarly, the terms in F2 can be estimated as 

11(1 _ i)(A„ + Vdivi;)|L < C||V2i;||2||a||^ < C6\\V^v\\2, 
p p 

||(i_i)(A^; + Vdivi;)||^i 
P P 

< - -){Av + VdWv)\\2 + ||V((i - -)iAv + Vdivw))||2 
P P P P 

< C(||a||oo||V2T;||2 + ||a|U||V3«||2 + ||Va|U||V2t;||2) 



<CS\\Vv\\h2, 



1 ^p'ia + p) p'ip) , 2{m + k) 2fc^^^||^ 



"fX a + p p 3{a + p) 3p 

7(||Va| 

< C,5||Va||2 



< C(||Va||2||a||^ + ||Va||2||m||^ 

2—p 2—p 
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,p'{a + p) 


P'ip) 


2(m 


+ fc) 


2k 


' a + /5 


P 


^ 3(a 


+ P) 




,p'{a + p) 


P'ip) 


2(m 


+ fc) 


2k. 


' a + p 


P 


^ 3(a 


+ P) 


~W 


p'ia + p) 


P'ip) ^ 


2(m 


+ k) 


2k. 


a + p 


P 


3(a 


^P) 


^Yp' 



<^(ll( 

+ liv[( 



12; 

< C(||aVa + mVa\\2 + ||aV^a + mV^a + (Va)^ + VaVm||2) 
<C5\\ya\\H2, 



^Vm||p < ||Vm||2||m + fc||6||^^||^ < C<5||Vm||2, 
c>7A 171 + k 3-2p 



II - 3^V™IIhi < ^^(l|Vm|l2 + llV^mlU) 

<C(||Vm||3||m + fc||6||^^|U + ||V2,7l||3||™ + ^||6|l^Tlloo) 

<CS\\ym\\H^. 
Thus we have 

\m,F2)\\LPnH^ <Cd\\{^a,Vv,Vm)\\H2. (3.3) 

Next, we estimate for || (F3, F4, i^5)||p. Ahuost as same as the estimation of || (Fi, F2)||/^pn/fi , 
we get 

||(F3, F4, i^5)||p < CSWiVa, Vf, V/i, Vm, Ve)||ff2 . (3.4) 
Inserting (|3.3p and (|3.4p into p.2p , we complete the proof of Lemma 13.11 □ 

Lemma 3.2. Let W ~ (U, h, m, e) be the solution to the problem (|2.ip . i/ien under 
the assumptions of Theorem if S > is sufficiently small then it holds that 

''^'^^^^ + llV^all^i + \\V\v,h,m,e)\\l2 < CS\\V{a,v, h,m, e)\\l, (3.5) 



dt 

where AI{t) is equivalent to || V(a, u, /i, to, e) ||^2; i.e., there exists a positive constant 
C2 such that 

C^^\\y{a,v,h,m,e)\\]j2 < M{t) < C2 |I V(a, w, /i, to, e) ||^2. (3.6) 

Proof. Let a be any multi-index with 1 < |a| < 3. Applying the operator 9" to 
(|2.ip and then taking inner product with 9"a, d'^v, d"h, d"m and d^e, one gets 

i^||a^"(a,„,/,,^,£)l|2 + A||va,"HI^ + A||div9>||2 + A||va,"(/i,TO,£)||2 

+ {d^e,d^F,) 

=■■ J? it) + j^t) + jnt) + Jtit) + j^it), 

where Jf^it), i = 1,2,3,4,5, are the corresponding terms in the above equation 
which will be estimated as follows. 
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Now let's estimate for J"(i). It follows from Lemma [2.31 that 

JI ^ {djdi{av^), djo) ^ {djdiav'' + diadjV^ + djadiv'' + didjv'^a, dja) 

< CIlV^allsllVallsMloo + il Va||2|| V«||oo + || V^^HsH Va||2||a|U 

< ^llVail^ + ^||V2a||2||^;||^, + ^||V2z;||^||a||l,., 



< -^llV^all^ + ydlV^all^ll^;!!^ + ||Va||^||Vz;||L + W^MlMl) 

where a = /? + ^, with |/3| = |^| = 1. 

For Jf, we have from multi-indices /3 and ^ with |/?| = 2, |^| = 1, 

Ji^^(9f+«9,(at;'),af+«a) 

< S\\V'a\\l + ydlV^all^llV^^;!!^ + + || V^z;|l2ll«llL) 

< ^llV^all^ + ji\\S7^a\\l\\v\\],. + \\Va\\l\\S7v\\l. + \\VMIMU- 



Hence 



J?<CS J2 \\9>\\l + CS J2 \\9>\\l (3-7) 

1<|q|<3 l<|a|<4 



75— ^Vm], i9jW) 



Similarly, by Lemma 12.31 and Holder inequality, can be estimated as 

,1 ,or.l Iw, ^ 1 t P'ja + p) P'jp) , 2(m + fc) 2fc 

■J2 - W (- - + Vdivw -( — — + - ^jVa 

p p a + p p i{a + p) 3p' 

2 
3"7A 

^ {dja{Av + Vdivu) + a{djAv + dj\/divv} + dja ■ Va + djui ■ Va + adj^a 
+ mdjVa + dj\7m, djv) 

< C(||Va||3||Vz;||6||V2«||2 + ||V3i;||2|lVi>!|2||a|U + || V(a, m)||3|l VaUeH Vt;||2 
+ l|V2a||2||Vt;||2||m|U + ||V2a||2||Vi>||2||a||oo + |lV2m||2|lVz;||3||Vm||2) 

< SiW^Ml + \\V'v\\l + l|V^(a,m)||2) + |(||Vz;l|2||(a,m)||^. + || Vz;||^. || Vm||^) 
<C6\\V\a,m)\\l + CS ^ \\dMl 

l<\a\<3 
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p p ^X^ a + p p 3[a + p) 3p' 

^ {Va{Av + Vdivw) + aV^v + VaVa + VaV77i + aV^a + mV^a + V^m, V^w) 

< C(||Va||3||V2i;||6||V3w||2 + llV^wll^llalloo + || Va||3|| Va||6|| V^t-lh 

+ ||Va||6||V\>||2||Vm||3 + ||V2a||2||V\>||2||(a,7n)|U + ||V2m||3||V\||2||m + fclle) 
<C6 J2 \\9>\\l + C5 J2 \\dS{v,m)\\l with a ^13 + \l3\^\i\^l, 

l<|a|<3 1<|q|<3 

Jl - (5f [(i - -)iAv + Vdiv.) - 1 (P>±^ + ^p±^ _ ^) Va 

^ ^ "^^V P -fX^ a + p p 3{a + p) 'ip' 

(aV^w + VaV^w + V^aV^w + aV^a + VaV^a + VmV^a + V'^mVa 

< CS\\V\\\l + f (||Va||L||V\-||^ + llV^^ll^llV^all^ + || V^a!|^||a||L 

+ \\V'a\\l\\V{a,m)\\l + || V^mH^H Va||L + II V=^a||^||m||L + || V^mbU Vm||2 
<CS ||a>||^ + C5 m{v,rii)\\l witha = /3 + t m^2, \^\ = l. 

1<|q|<3 l<|a|<4 

Incorporating the above estimates, it holds that 

Jnt)<CS Y I|5"«ll2 + C^ Y Il^"(^''"^)ll2- (3-8) 

1<|q|<3 l<|a|<4 

Moreover, J^{t), J^{t) and J^{t) can be estimated similarly by using Lemma 
12.31 and Holder inequality, 

\J^{t),J2{t),J^{t)\<C5 Y I|5:«ll2 + ^^^ E \\d:{v,h,m,e)\\l (3.9) 

1<|q|<3 l<|a|<4 

On the other hand, we apply to (|2.ip „ with 1 < |q;| < 2, and then take inner 
product with V9"a to yield 

{dyuydSa)+^\\Vdy\\l = X{dSAv + dSVd\vv,dSVa) + {d^F^.d^Va), (3.10) 

and similarly from (|2.ip -^. 

{d^v, Vd^at) = -7(5^«, Va,"div«) + {dy, Vd^F,). (3.11) 

Adding p.lOp and p. lip together implies 

+ {dSF,,asva) - -/{oy, va^divv) + {ay, va^F,). 
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By Holder inequality and similar to the estimation of J2it), the right hand side 
can be bounded by 

|||V5>||^ + C J2 \\d:^v\\l.+CS J2 \\d>\\l + CS J2 \\d^iv,h,m,s)\\l 

1<|q|<2 l<|a|<3 1<|q|<4 

(3.13) 

Therefore, if we define 

M{t) = C, \\d:ia,v,h,m,s)\\l+ ^ (a>,V9,%), 

1<|q|<3 l<|a|<2 

and choosing 6 sufficiently small, then p.7p - p.9p and (|3.13p imply that 
dM{t) 



dt 



+ Ci{\\y'a\\jj, + \\V'{v,h,m,e)\\j,,) < CS\\y{a,v,h,m,s)\\i 



where Ci is a positive constant independent of 6. Thus we arrive at the proof of 
the lemma. □ 

4. Optimal convergence rates 

The optimal convergence rates can be proved by first improving the estimates 
given in Lemma |3. II and Lemma 13.21 to the estimates on the L^-norms of solutions 
to higher power and then letting the power tend to infinity. By the inequalities 
p.ip and (|3.5p . we have the following lemmas. 

Lemma 4.1. Let W = {U,h,m,e) be the solution to the problem (|2.ip . then under 
the assumptions of Theorem if S > is sufficiently small, then for any integer 
n > 1, and for some p G f )■ holds that 

t pt 

^ II WT^"// II 2n j„ ^ / TP \2n i //^r\2n / /i i „\^IIv72t 



(l + sy||VW^(s)||f rfs< (C£;o)'" + (C(5)^" / (l + s)'||V^W^(s)||^"i&, (4.1) 
Jo 

where I ~ 0,1, ■ ■ ■ , N ~ [2n{j^ ~ ~ 2], the constant Eq is given in Lemma \S.l\ 

Lemma 4.2. Let W = {U,h,m,e) be the solution to the problem (|2.ip . then under 
the assumptions of Theorem \1.2l if S > is sufficiently small, then for any integer 
n > 1, and for some p e [1; f ) ■ holds that 

(1 + tyM{tr + n f\i + syM{sr-'\\v^wis)\\i, ds 



<2M(0)" + (C3£;o)'" + 8C2n(|--i) / {1 + sy~^M{sr~^\\V'Wis)Wj,,ds, 

(4.2) 

where I = 0, 1, • • • ,N= [2n(^ — j) — 2], the constant C2 is given in Lemma \3.'A 
C3 is independent of 5. 

Remark 4.3. Lemma |4?T] and Lemma |4?2] arc similar to that in [12]. But Lemma |4T 



and Lemma [4.21 in this work hold for general p with p S [1, f )• Note that lemmas 
hold only for p = 1 in [T^]. For completeness, we state the proofs of Lemma [4.11 
and Lemma[4.2las follows. 
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Proof of Lemma \4-l\ Fix any integer n > 1. By taking p.ip to power 2n and 
multiplying it by (1 + t)', / = 0, 1, ■ • • , N, integrating the resulting inequality over 
[0, t] gives that 



(4.3) 



/ (l + T)'||VVF(T)||2"dT < (C£;o)'" / (l + T)-2»(|;-3)+'dr 
Jo Jo 

+ iCSf" f{l + T)'[ r(l + ^-'S)-(*-^)|lVW^(s)||ff2ds]'"dr. 
Jo Jo 

It follows from the Holder inequality that 

(l+T-s)-(*-3)||VW^(s)||^2ds]'" 

<[( (l + T-s)-''i(l + .s)-''Ms]^"~^ (4.4) 
Jo 

X [\l + T- + s)'||VM^(s)||^Ws, 



where 



and 



,3 1 2 , 2n 

ri = ( 



2p 4 3n'2n-l 

f2 



2n - 1 

Notice that | - ^ < ri < |j - i and ra € [0, ri] for n > 1 and < Z < TV = 
[2n(^ ^ i) ^ 2], from Lemma one deduces that 

(1 + r - .s)^'-! (1 + s)-''Ms < Ci(ri, r2)(l + r)"''^ < C(l + T)-'■^ (4.5) 



where Ci(ri,r2) given by (|2.8I) is bounded uniformly for n > 1. Hence, (|4.3p 
together with (|4.4p and (j4.5p leads to 

r(l + r)'||VW^(r)||rdr<(Cii;o)^"— 3 \ 

Jo 2n(^ - ^) - Z - 1 

+ iCdf'' J\l + sy\\\/W{s)ffi\ J\l + T- s)-^drds 
< (CEo)^'' + (CSf'' f (l + .s)'||VVF(s)||2,",ds 



(4.6) 



< (C-Bo)'" + (C(5)2" / {l + sy{\\VW{s)\\l" + \\V^W{s)\\j^^)ds. 
Jo 

Here we have used the fact 

2n(|;-i)-i-l>2n(|.-i)-2,.(|.-i) + 2-l.l, 

Thus if (5 > is sufficiently small such that (C7i5)^" < ^ in the final inequality of 
(HH), then (gH) implies (gll]). We finish the proof of Lemma O □ 
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Proof of Lemmall^ Multiplying ([33]) by n{l + tyM{t)"-^ ioi I = 0,1, ■ ■ ■ , N 
and integrating it over [0,t] give that 

(1 + t)'A/(t)" + n f\l + s)'A/(s)"-i|l V2W'(s)||^ids 
Jo 

<MiOy' + CSn [ (l + s)'M(s)"-i||VTy(s)||^ds + Z [ (1 + s)'"^M(s)"ds. 

(4.7) 

For the second term on the right hand side of (j4.7p . from the Young inequality, 
(|3.6p and Lemma WA\ it holds that for any rj > 0, 



Sn /*(l + -5)'^^(s)""'llVW(s)|l 
Jo 



?ds 



< 6n f\l + sYi^^Misr + -^\\VW{s)\\l-]ds 
Jo n nri'^ ^ 

t 

n-li 



<SnC2V / (l + s)'M(s)"-^(||VM^(s)||^ + ||V^M/(.s)||^Od^ 
Jo 

+ [\i + sy\\ww{s)\\f^ds 

Jo 

<SnC27] [ (l + s)'Af(s)"-i||VVF(s)||^ds 
Jo 

+ 5nC2V I (l + s)'M(s)"-'||V2M^(s)||^ids 
Jo 

+ ,5771-" [(CSo)'" + j\l + syW^^W{s)\\%ds\ 



(4.8) 



Jo 

<Srj^~''{CEo)^" + SnC27l [ (1 + s)'M(s)"-i||Viy(s)||2ds 

Jo 

+ SnC2V I {l + sfM{sY~^\\V^W{s)\\]J^ds 
Jo 

+ 5i]^~'\C5f'^Gi~^ f (l + s)'Af(s)"-i||v2Ty(s)|l^ids. 
Jo 



Choosing rj = in (|4.8p . it holds that 



Sn (l + s)'A/(s)"-^||VM^(s)||^ds 





< 2Si2C2y'-\CEaf'' + Sn[l + -{CSf''{2Cl)''-'] f\l + sY M{sr-^\\V^W{s)\\jj,ds 

^ Jo 

\2ji I x„^^ I t/^X\'2n} / ^i , „\l A/rt „\n-l\\T-,2j 



<d{CEor' +6n[l + {C5y''] / (1 + s)'A^(s)"-Ml V'VK(s)||^ids. 

(4.9) 
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Thus if (5 > is sufficiently small such that C5 < 1 in (gH), then (CfJ)^" < 1 for 
any n > 1. And from (|4.9|) . it is easy to show that 



^ (4.10) 

<{CEof'' + 2Sn [ {1 + syM{sr-^\\V^W{s)\\l^ds. 



Similarly, we can estimate for the third term on the right hand side of (14.71) as 

I f (l + s)'-iAf(s)"ds 
Jo 

<[CEoY'' + 5n{C5Y'^^^ f (l + s)'M(s)"-i||V2M^(s)||^ids (4.11) 

Jo 

+ 2IC2 I {l + sf-^M{sy'-^\\V^W{s)fH.ds, 
Jo 

which together with (|4.7|) and ()4.10p arrives at 

(1 + t^MitY + n f\l + s)'M(s)"-i II v2VF(s) lipids 



< Af(0)" + (C£;o)'" + (5n[C+(C(5)""-^] / (1 + s)'Af (s)"-^||V"t^(.s)|||^ids 

+ 2IC2 f {i + sy-'M{sy'-'\\v^w{s)\\l^ds. 

Jo 

(4.12) 

Choose 6 > sufficiently small such that for any 71 > 1, it follows that 

5[C + C,5)2"-i] < i, 
then it follows from (|4.12p that 

(1 + t)'M(i)" + n f\l + s)'M(s)"-i||V2W(s)||^ids 
Jo 

< 2M(0)" + (C3S0)'" + 4ZC2 / {1 + sy-HI{sr~^\\V^W{s)\\jJ^ds, 



which gives ()4.2p because I < N < 2n{^ ~ j). Thus we complete the proof of 
Lemma 14.21 □ 



Proof of Theorem \1.2\ Let (5 > be small enough such that Lemma |4?2] holds for 
any n>2. For any fixed integer n > 2, from Lemma 14.21 '^s get that the inequality 
(|4^ holds for any ^ = 0, 1, • • • , ^. When I = 1, (|42]) reads 



t 

n-l II v72ti^/„m|2 



(1 + t)A/(t)" + n / {l + s)M{sT''\\V'W{s)\\ij,As 

^ (4.13) 

< 2A/(0)" + (C3S0)'" + 8C2n(|- - i) / Af(s)"-i||V2w^(s)||2,,ds. 

2p 4 Jo 



(4.14) 
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It follows from ()1.3p in Proposition 11.11 that 

f M{sr-'\\V^W{s)\\l^ds < [snpM{s)r-' T ||V'W^(s) Ill/ids 

JO s>0 Jo 

which together with (|4.13p shows that 

{l+t)M{t)" +n f {1 + s)M{sy'-^\\V^W{s)\\l^ds 
Jo 

< 2M{0r + {C^Eof^' + 8C2n(|- - ^){C2CoS^r. 

2p 4 

For 1 < ^ < A'^, by induction one can arrive at 

(1 + t)'Af(i)" + n f\l + syM{sr-'\\V^Wis)\\l^ds 
Jo 

< [2M(0)" + iCsEof-] jySC^i^ + n[8C2i^ ~ l^^C^CoS'T. 

In fact, suppose that (|iTT5)) holds for 1 < / < - 1. Then from it holds that 

(1 + ty+Hi{ty' + n f\l + sy+Hi{sy'-^\\v^w{s)\\l^ds 

Jo 

< 2Af (0)" + (C3£;o)'" + 8C2n( A - i) / (1 + sy M{sr~'\\\7^W{s)\\j,^ds 

< [2A/(0)" + {C^Eof-] + 8C2(|- - i){[2M(0)" + (CgiJo)'"] Et^^^lf - j)]-' 
+ n[8C2{^-l)nC2CoS'r} 



(4.15) 



(4.16) 



which combining with (|4.14p gives that (|4.15p holds for any 1 < Z < iV. 
Specially, 



,[ 8g2(§-i)r-i 3 1. 



(1 + tfM{tr < [2Af(0)" + (C3go)'"] ^^ Ts i\ . + «[8C'2(7e - i)f{C2Co5^r 



Note that 



M^-J)-3<^=[2n(|-4)-2]<2.(l-i)-l. 



It is not difficult to prove that 

which gives that 

A/(i)^ < C[A^(0)" + £;^" + 5^"]^(l + i)"^*"^^^^- 
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Since M(0), Eq and 6 are independent of n, one gets 

[A'/(0)" + E^" + S"^"]^ max{VAf(0), Eq, 6}, as n ^ oo. 
The above relation implies that 

M(t)^ <Cma.x{^/mO),Eo,d}il+ty^^~i\ 

that is, 

<Cmax{yM(0),i?o,<5}(l + t)-'^(f'2;i)^ 

which together with Lemma 12.31 implies (|1.7p and p.8p . 

Now, estimate for p.6p . For this purpose, applying (|1.8p . Lemma [^TT] and Lemma 
O leads to 

\\W{t)h < CEoil + t)--(f'2;0) f\l + t~ s)--'^P''-''\\\F{W)\\p + ||^^(W^)||2)ds 

Jo 

< CEoil + ty'-P'^''^^ +CS [ (l+t-s)-"(P^2;0)||yp^(^)||^^j^ 

Jo 
Jo 

< C(l 4-t)-'"(P^^'°^ 

By interpolation, we have that (|1.6p holds for any 2 < 5 < 6. 
For (jOj) . from ([2TT|) . we get 

\\dtW{t)\\2 < WlAiYvh + ^^2, ^^3,-^4,-^5)112 + II - 7Va + \Av + AVdivi;||2 
+ ||AA/i||||2 + ||AAm||2 + ||AAe||2 

< C(||V(a,w,/i,m,£)||2 + llV2(u,/i,TO,£)l|2) 

< CEoil +ty'^P''^'^\ 

Thus, (|1.9p is proved. The proof of Theorem 11.21 is complete. □ 
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